Modeling Neutrino Quasielastic Cross Sections 
on Nucleons and Nuclei 



A. Bodek*, H. Budd* and J. Arrington^ 



m 
O 
O 
(N 

Oh! 

(N 
> 

O 

a^ 
o 

m 

o 

D 
D 



X 



* University of Rochester, Rochester, NY 
^Argonne National Laboratory, Argonne, Illinois 60439, USA 



Abstract. We calculate the total and differential quasielastic cross sections for neutrino and an- 
tineutrino scattering on nucleons using up to date fits to the nucleon elastic electromagnetic form 
factors G^, G'^, G^, G^, and Fa and pseudoscalar form factors. We compare predictions of the cross 
sections for nucelons and nuclei to experimental data. (Presented by Arie Bodek at CIPANP2003, 
New York City, NY 2003) 



1. INTRODUCTION 

Experimental evidence for oscillations among the three neutrino generations has been 
recently reported [1]. Since quasielastic (QE) scattering forms an important component 
of neutrino scattering at low energies, we have undertaken to investigate QE neutrino 
scattering using the latest information on nucleon form factors. 

Recent experiments at SLAC and Jefferson Lab (JLab) have given precise measure- 
ments of the vector electromagnetic form factors for the proton and neutron. These form 
factors can be related to the form factors for QE neutrino scattering by conserved vec- 
tor current hypothesis, CVC. These more recent form factors can be used to give better 
predictions for QE neutrino scattering. 

The hadronic current for QE neutrino scattering is given by [2] 



< p{p2)\Jt\n{pi) >-- 
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where q = ky — k^,, t, = [jip — \) — jin, and M = {nip -\- m„)/2. Here, jip and jln are 
the proton and neutron magnetic moments. We assume that there are no second class 
currents, so the scalar form factor Fy and the tensor form factor need not be included 
Using the above current, the cross section is 
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FIGURE 1. The QE neutrino cross section along with data from various experiments. The calculation 
uses Ma=1.00 GeV, gA=~1.267, M^=0.71 GeV^ and BBA-2003 Form Factors. The soHd curve uses no 
nuclear correction, while the dashed curve [5] uses a Fermi gas model for carbon with a 25 MeV binding 
energy and 220 Fermi momentum. The dotted curve is the prediction for Carbon including both Fermi gas 
Pauli blocking and the effect of nuclear binding on the nucleon form factors [8]. The data shown [3] are 
from FNAL 1983, ANL 1977, BNL 1981, ANL 1973, SKAT 1990, GGM 1979, LSND 2002, Serpukov 
1985, and GGM 1977. 
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Although we have have not shown terms of order (nii/M)^, and terms including 
Fp{q^) (which is multiplied by (mi/M)^), these terms are included in our calcula- 
tions [2].) The form factors Fy{q^) and t,Fy{q^) are given by: 



1 



4M^ 



1 



4M^ 



We use the CVC to determine G\{q-^) and G\j{q^) from the electron scattering form 
factors G^Eiq^), Gl{q^), G^jiq^), and G^(^2). 



GW) = Gl{q^)-Gl{q% G\{q^) = G^^q^) - G],{q\ 
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The axial form factor Fa and the pseudoscalar form factor Fp (related to Fa by PCAC) 
are given by 

Ml) 

In the expression for the cross section, Fp{q^) is multiplied by {m[/M)^. Therefore, in 
muon neutrino interactions, this effect is very small except at very low energy, below 
0.2 GeV. FA{q^) needs to be extracted from QE neutrino scattering. At low Q^, FA{q^) 
can also be extracted from pion electroproduction data. 

Previously, people have assumed that the vector form factors are described by the 
dipole approximation. 



GD{q^) = - ^— M^ = 0.nGeV^ 



G^E = GD{q'), G| = 0, G^M = HpGD{q'), G^, = MnGz)(?"). 

We refer to the above combination of form factors as 'Dipole Form Factors'. It is an 
approximation that has been improved by us in a previous publication [3]. We use our 
updated form factors to which we refer as 'BBA-2003 Form Factors' (Budd, Bodek, 
Arrington). We also use our updated value [3] of Ma 1.00 ± 0.020 GeV which is in 
good agreement with the theoretically corrected value from pion electroproduction [4] 
of 1.014 ±0.016 GeV. 



2. COMPARISON TO EXPERIMETAL DATA 

Figures 1, 2, and 3 show the QE cross section for v and V with BBA-2003 Form Factors 
and Ma=1.00 GeV. The normalization uncertainty in the data is approximately 10%. The 
solid curve uses no nuclear correction, while the dashed curve [5] uses a NUANCE [6] 
calculation of a Smith and Moniz [7] based Fermi gas model for carbon. This nuclear 
model includes Pauli blocking (see Figure 4(a) ) and Fermi motion, but not final state 
interactions. The Fermi gas model was run with a 25 MeV binding energy and 220 MeV 
Fermi momentum. The dotted curve is the prediction for Carbon including both Fermi 
gas Pauli blocking, and the effect of nuclear binding on the nucleon form factors as 
modeled by Tsushima et al [8] (see Figure 4(b)). Note that this model is only valid for 
less than 1 GeV^, and that the binding effects on the form factors are expected to 
be very small at higher Q^. Both the Pauli blocking and the nuclear modifications to 
bound nucleon form factors reduce the cross section relative to the cross section with 
free nucleons. 

The updated form factors improve the agreement with neutrino QE cross section data 
and give a reasonable description of the cross sections from deuterium. 

We plan to continue to study the nuclear corrections, adopting models which have 
been used in precision electron scattering measurements from nuclei at SLAC and 
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FIGURE 2. Same as Figure 1 with the Ey axis limit changed to 2 GeV. 



JLab. For example, we plan to study the Pauli blocking correction using an improved 
Fermi Gas model with a high momentum tail [9], as well as more sophisticated nuclear 
spectral functions. In addition, we will continue to update the extraction of Ma from 
previous neutrino experiments, using the updated versions of the input parameters and 
electromagnetic form factors. 

This work is supported in part by the U. S. Department of Energy, Nuclear Physics Di- 
vision, under contract W-31-109-ENG-38 (Argonne) and High Energy Physics Division 
under grant DE-FG02-91ER40685 (Rochester). 
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FIGURE 3. The QE antineutrino cross section along with data from various experiments. The calcu- 
lation uses Ma=1.00 GeV, gA=-l.267, M^=0.71 GeV^ and BBA-2003 Form Factors. The soHd curve 
uses no nuclear correction, while the dashed curve [5] uses a Fermi gas model for carbon with a 25 MeV 
binding energy and 220 MeV Fermi momentum. The dotted curve is the prediction for Carbon including 
both Fermi gas Pauli blocking and the effect of nuclear binding on the nucleon form factors [8]. The data 
shown are from SKAT 1990, GGM 1979, Serpukov 1985, and GGM 1977. 
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FIGURE 4. (a) The Pauli blocking suppression for a Fermi gas model for carbon with a 25 MeV 
binding energy and 220 MeV Fermi momentum, (b) The ratio of bound to free nucleon form factors for 
Fi, F2, and fromref [8]. Note that this model is only valid for less than 1 GeV^, and that the binding 
effects on the form factors are expected to be very small at higher Q^. 
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